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Elliott H. Lieb and Robert Seiringer 
Department of Physics, Jadwin Hall, Princeton University, P. O. Box 708, Princeton, New Jersey 08544 



vn ■ 
o ■ 
o : 

(N . 

X) : 

(N : 

> • 
m ■ 
(N ■ 
O 
(N 

o . 



X 



Jakob Yngvason 

Institut fiir Theoretische Physik, Universitdt Wien, Boltzmanngasse 5, A-1090 Vienna, Austria 

(Dated: Feb. 28, 2005) 

The validity of substituting a c-number z for the k = mode operator Oq is established rigorously 
in full generality, thereby verifying one aspect of Bogoliubov's 1947 theory. This substitution not only 
yields the correct value of thermodynamic quantities like the pressure or ground state energy, but 
also the value of \z\'^ that majcimizes the partition function equals the true amount of condensation 
in the presence of a gauge-symmetry breaking term — a point that had previously been elusive. 
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One of the key developments in the theory of the Bose 
gas, especially the theory of the low density gases cur- 
rently at the forefront of experiment, is Bogoliubov's 
1947 analysis jlj of the many-body Hamiltonian by means 
of a c-number substitution for the most relevant opera- 
tors in the problem, the zero-momentum mode operators, 
namely *. Naturally, the appropriate 

value of z has to be determined by some sort of consis- 
tency or variational principle, which might be compli- 
cated, but the concern, expressed by many authors over 
the years, is whether this sort of substitution is legit- 
imate, i.e., error free. We address this latter problem 
here and show, by a simple but rigorous analysis, that it 
is so under very general circumstances. 

The rigorous justification for this substitution, as far 
as calculating the pressure is concerned, was done in a 
classic paper of Ginibre 0] in 1968, but it does not seem 
to have percolated into the general theory community. 
In textbooks it is often said, for instance, that it is tied 
to the imputed 'fact' that the expectation value of the 
number operator no = a^an is of order V — volume. 
(This was the argument in 1]). That is, Bose-Einstein 
condensation (BEG) justifies the substitution. As Gini- 
bre pointed out, however, BEG has nothing to do with 
it. The z substitution still gives the right answer even 
if no is small (but it is a useful calculational tool only if 
no is macroscopic). Thus, despite and the thorough 
review of these matters in 0, there is some confusion in 
the literature and clarification could be useful. 

In this short note we do three things. 1.) We show 
how Ginibre's result can be easily obtained in a few sim- 
ple lines. While he used coherent states, he did not use 
the Berezin-Lieb inequality 0, H, j derived later, which 
efficiently gives upper bounds. This inequality gives ex- 
plicit error bounds which, typically, are only order one 
compared to the total free energy or pressure times vol- 
ume, which are order N — particle number. 

2.) This allows us to go beyond 2] and make c-number 
substitutions for many k-modes at once, provided the 
number of modes is lower order than N . 



3.) We show how the optimum value of z yields, in 
fact, the expectation value (no) in the true state when a 
gauge breaking term is added to the Hamiltonian. More 
precisely, in the thermodynamic limit (TL) the |zp that 
maximizes the partition function equals |(ao)P and this 
equals (rio), which is the amount of condensation — a 
point that was not addressed in full generality in previ- 
ous work 0, 0j 0- The second of these equalities has 
previously only been treated under some additional as- 
sumptions or for some simplified models 0, 1^ ■ 

While we work here at positive temperature k^T = 
1//3, our methods also work for the ground state (and 
are even simpler in that case). To keep this note short 
and, hopefully, readable, we will be a bit sketchy in places 
but there is no difficulty filling in the details. 

The use of coherent states 0, to give accurate up- 
per and lower bounds to energies, and thence to expecta- 
tion values, is effective in a wide variety of problems [l^ . 
e.g., quantum spin systems in the large S limit [a, the 
Dicke model [l^, the strong coupling polaron [lj|, and 
the proof that Thomas- Fermi theory is exact in the large 
atom limit 0, 0| . For concreteness and relevance, we 
concentrate on the Bose gas problem here, and we dis- 
cuss only the total, correct Hamiltonian. Nevertheless, 
the same conclusions hold also for variants, such as Bo- 
goliubov's truncated Hamiltonian (the "weakly imperfect 
Bose gas" 0,0) or other modifications, provided we are 
in the stability regime (i.e., the regime in which the mod- 
els make sense) . We are not claiming that any particular 
approximation is valid. That is a completely different 
story that has to be decided independently. The method 
can also be modified to incorporate inhomogeneous sys- 
tems. The message is the same in all cases, namely that 
the z substitution causes no errors (in the TL), even if 
there is no BEG, whenever it is applied to physically sta- 
ble systems. Gonversely, if the system is stable after the 
z substitution then so is the original one. 

We start with the well-known Hamiltonian for bosons 
in a large box of volume V, expressed in terms of 
the second-quantized creation and annihilation operators 
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^k' '^k satisfying the canonical commutation relations, 
H = Y^ fc^flkflk + ^Y1 ^(P)«k+paq-pakaq: (1) 



2V 

k k.p.q 

(with h — 2m = 1). Here, v is the Fourier transform of 
the two-body potential v{r). We assume that there is a 
bound on the Fourier coefficients |j^(k)| < (p < oo. 

The case of hard core potentials can be taken care of 
in the following way. First cut off the hard core potential 
V a.t a height 10^^ eV. It is easy to prove, by standard 
methods, that this cutoff will have a negligible effect on 
the exact answer. After the cutoff ip will be about lO'^^ 
eV , and according to what we prove below, this sub- 
stitution will affect the chemical potential only by about 
(p/V, which is truly negligible when V = 10^'^ A"^. 

If we replace the operator Oq by a complex number z 
and Oq by z* everywhere in H we obtain a Hamiltonian 
H'{z) that acts on the Fock-space of all the modes other 
than the mode. Unfortunately, H'{z) does not com- 
mute with the particle number iV> = X]k/o^k"k- 
convenient, therefore, to work in the grand canonical en- 
semble and consider i/^ = H — fiN = H — ^{aQaQ + N^) 
and, correspondingly, H'^{z) — H'{z) — -I- iV>). 

The partition functions are given by 

e/3Vp(M) ^ = Tr„ cM-W (2) 
^fSVp'M ^ ^ J T^^, exp[-/3i/;(z)] (3) 

where H is the full Hilbert (Fock) space, H' is the Fock 
space without the Qq mode, and d^z = dxdy with 
z — X + iy. The functions p{fJ.) and p'{fJ.) are the corre- 
sponding finite volume pressures. 

The pressure p{fJ.) has a finite TL for all fi < /icriticai, 
and it is a convex function of /i. For the non-interacting 
gas, /^critical = 0, but for any realistic system ^critical = 
-|-oo. In any case, we assume /j, < ^critical, in which case 
both the pressure and the density are finite. 

Let \z) = exp{— |zp/2 -I- zog} |0) be the coherent state 
vector in the Og Fock space and let 11(2) = \z){z\ be 
the projector onto this vector. There are six relevant 
operators containing Oq in H^^ which have the following 
expectation values 10] (called lower symbols) 



{z\ao\z) 
{z\al\z) 



[z\aoao\z) 
{z\aoao\z) 



(z|aJjao|z) ^ \zY 

* I 



Each also has an upper symbol, which is a function of 
z (call it u{z) generically) such that an operator F is 
represented as F ^ J d'^z u{z)Il{z) . These symbols are 



J,gUQ 



^ , "0"0 
z , agagflQaQ 



U|4-4|z|2 + 2. 



It will be noted that the operator H'^{z), defined above, 
is obtained from iJ^ by substituting the lower symbols 



for the six operators. If we substitute the upper symbols 
instead into we obtain a slightly different operator, 
which we write as H'^{z) = H'^{z) + 5fj,{z) with 



5,{z)=^i+^[{-A\z\^ + 2)v{G) 



-5]a;;ak(2KO) + Kk) + K-k))]. (4) 

ky^O 

The next step is to mention two inequalities, of which 
the first is 



S(m) > S'(m) 



(5) 



This is a consequence of the following two facts: The 
completeness property of coherent states, / d\\i{z) = 
Identity, and 

(6) 

where cj) is any normalized vector in TC . This is Jensen's 
inequality for the expectation value of a convex function 
(like the exponential function) of an operator. 

To prove lO we take ^ in © to be one of the nor- 
malized eigenvectors of H'^{z), in which case exp{((/)| — 
/3Hl{z)\<P)} = {(t>\e^p{-l3H'^{z)}\(t>). We then sum over 
all such eigenvectors (for a fixed z) and integrate over z. 
The left side is then S(/i), while the right side is S'(/i). 

The second inequality |j, S la| is 



(7) 



where S"(/i) is similar to S'(/i) except that i?^'i(z) is re- 
placed by H'^{z). Its proof is the following. Let e Ti 
denote the complete set of normalized eigenfunctions of 
H^. The partial inner product \^ j{z)) — {z\^j) is a vec- 
tor in Ti.' whose square norm, Cj{z) = j{z)\^ j{z))n' 
satisfies J d'^zcj{z) = 1. By using the upper symbols, 
we can write {^j\H^\<^j) = J d^z (■i'j{z)\H'/^{z)\-9j{z)) = 
/ dh {^'^{z)\Hf{z)\'^'^{z))c/{z), where |*;-(z)) is the nor- 
malized vector Cj{z)~'^^^'ij{z). To compute the trace, we 
can exponentiate this to write S(/^) as 

^expj-/? J d'zc,{z)(mz)\H';^iz)\mz))Y 

Using Jensen's inequality twice, once for functions and 
once for expectations as in JHl, S(/i) is less than 

j dhc,{z)e^v{{%{z)\~ PH;{z)\%{z))] 
<J2 [ d'zc,iz){mz)\e^p{-PH';iz)}\%iz)). 

Since Trn(z) = 1, the last expression can be rewritten 
j d\ ^($, |n(z) ® exp{-/3i7;(z)} 1$,) = S"(m). 
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Thus, we have that 

S'(^) < S(^) < (8) 

The next step is to try to relate to To this 

end we have to bound (5p(z) in This is easily done in 
terms of the total number operator whose lower symbol 
is iV'(z) = |zp + X^k^^o '^k'^k- terms of the bound (p 
on v{p) 

|^^(z)| <2^(7V'(z) + i)/F+H . (9) 

Consequently, and are related by the in- 

equality 

S"(/i) < E'in + 2^/V^)e''(l^l+^/^). (10) 

Equality of the pressures p(/i), p'(a*) and p" {^i) in the TL 
follows from ||SJ| and ((Tn|l . 

Closely related to this point is the question of relat- 
ing 2(/i) to the maximum value of the integrand in ((SJ, 
which is max^ Tr^' exp[-/3i?^(z)] = e^^^P""". This latter 
quantity is often used in discussions of the z substitution 
problem, e.g., in refs. One direction is not hard. 

It is the inequality (used in ref. Q) 

S(/^) > maxTr„, exp[-/3i/;(z)], (11) 

and the proof is the same as the proof of |SJ|, except that 
this time we replace the completeness relation for the 
coherent states by the simple inequality Identity > n(z) 
for any fixed number z. 

For the other direction, split the integral defining 
S"(/i) into a part where |zp < ^ and |zp > ^. Thus, 

S"(m) < CniaxTr^. exp[-/3i/;(z)] 

+ - I d2z|z|2Tr„,exp[-/3i7;(z)]. (12) 

Dropping the condition |zp > ^ in the last integral and 
using |zp < N'{z) = N"{z) + 1, we see that the second 
line in (O is bounded above by ^~^E" {fi)[V p" (n) + 1], 
where p"{n) denotes the density in the H'^ ensemble. 
Optimizing over ^ leads to 

E"ifi) < 2[Vp"{^i) + 1] maxTr„, exp[-pH';{z)]. (13) 

Note that p"{^i) is order one, since p"(/z) and p{p) agree 
in the TL (and are convex in p), and we assumed that 
the density in the original ensemble is finite. By JHl, 
H'li > H'^+2v/v - ImI - V/V. and it follows from 0, ^ 
and pi(l that p™^^ agrees with the true pressure p in 
the TL. Their difference, in fact, is at most OiluV /V). 
This is the result obtained by Ginibre in Q by more 
complicated arguments, under the assumption of super- 
stability of the interaction, and without the explicit error 
estimates obtained here. 



To summarize the situation so far, we have four ex- 
pressions for the grand canonical pressure. They are all 
equal in the TL limit, 

p{p)=p'{p)^p"{p)^p^^--{p) (14) 

when p is not a point at which the density can be infinite. 

Our second main point is that not only is the z sub- 
stitution valid for Oq but it can also be done for many 
modes simultaneously. As long as the number of modes 
treated in this way is much smaller than N the effect 
on the pressure will be negligible. Each such substitu- 
tion will result in an error in the chemical potential that 
is order tp/V . The proof of this fact just imitates what 
was done above for one mode. Translation invariance is 
not important here; one can replace any mode such as 
X^k 5k'^k t>y a c-number, which can be useful for inhomo- 
geneous systems. 

A more delicate point is our third one, and it requires, 
first, a discussion of the meaning of 'condensate fraction' 
that goes beyond what is usually mentioned in textbooks, 
but which was brought out in P, 0, 0|. The 'natu- 
ral' idea would be to consider V~^{nQ). This, however, 
need not be a reliable measure of the condensate frac- 
tion for the following reason. If we expand exp{— 
in eigenfunctions of the number operator no we would 
have (no) = Sn^7(^)i where 7(n) is the probability 
that no = n. One would like to think that 7(n) is sharply 
peaked at some maximum n value, but we do not know 
if this is the case. 7(n) could be flat, up to the maximum 
value or, worse, it could have a maximum at n = 0. Re- 
call that precisely this happens for the Heisenberg quan- 
tum fcrromagnet |l7j : by virtue of conservation of total 
spin angular momentum, the distribution of values of the 
z-component of the total spin, S*^, is a strictly decreasing 
function of IS*^]. Even if it were flat, the expected value 
of would be half of the spontaneous magnetization 
that one gets by applying a weak magnetic field. 

With this example in mind, we see that the only phys- 
ically reliable quantity is \\m\^Q\\mv^ooV~^ {no) ^^\, 
where the expectation is now with respect to a Hamil- 
tonian H^,,x = H^, + ^/V{\aQ + A*a5) 1]. Without 
loss of generality, we assume A to be real. We will 
show that for almost every A, the density ^{V po) con- 
verges in the TL to a (5-function at the point = 
limy_>oo l^maxP/^, where Zmax maximizes the partition 
function Tr-^' exp{— \{'^T\- That is, 

^"'("0),,.A = l^~'|(ao)M,A|' = l^^'kmaxP (15) 

in the TL. This holds for those A where the pres- 
sure in the TL is differentiable; since p{p^ A) is convex 
(upwards) in A this is true almost everywhere. The 
right and left derivatives exist for every A and hence 
limA-^o-H limy^oo y~^\{ao) exists. 

The expectation values {no)fi,\ and (a(,)fi,\ are ob- 
tained by integrating (|zp — 1) and z, respectively, with 
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the weight W^,xiz) = E{n,X)-^TTn'{z\eM-PH^,x}\z). 
We will show that this weight converges to a ^-function 
at ^max in the TL, implying p5|l . If we could replace 
Wf,,x{z) by W^,o(z)e-''^^(^+^*', this would follow from 
Griffiths' argument [13 (see also [H Sect. 1]). Be- 
cause [i?, Oq] ^ 0, W^,A is not of this product form. 
However, the weight for A), which is W'^ ^{z) = 

E"{fi,\)'^TTn' exp{-PH'^ -^{z)}, does have the right 
form. In the following we shall show that the two weights 
are equal apart from negligible errors. 

Equality ((njl holds also for all A, i.e., p{fJ.,X) — 
p"(^,A) — p™^^{^,\) in the TL. In fact, since the up- 
per and lower symbols agree for Oq and Oq, the error 
estimates above remain unchanged. (Note that since 
VV\aa + a*o\ < S{N + ^) + V/S for any 5 > 0, p{n, A) 
is finite for all A if it is finite for A = in a small interval 
around /i.) At any point of differentiability with respect 
to A, Griffiths' theorem 17] (see Igj, Cor. 1.1]), applied to 
the partition function S"(A, /x), implies that ^^{CVV) 

converges to a (5-function at some point C on the real axis 
as 1^ ^ oo. (The original Griffiths argument can easily 
be extended to two variables, as we have here. Because 
of radial symmetry, the derivative of the pressure with 
respect to ImA is zero at any non-zero real A.) More- 
over, by comparing the derivatives of p" and p"^^^ we see 
that C = limy^oo z^ax/VV, since ^max/v^ is contained 
in the interval between the left and right derivatives of 
piiiaxj^^^ A) with respect to A. 

We shall now show that the same is true for W^^.a- 
To this end, we add another term to the Hamiltonian, 
namely eF = eV J d'^zn{z)f{zV~^/'^), with e and / 
real. If /(C) is a nice function of two real variables with 
bounded second derivatives, it is then easy to see that 
the upper and lower symbols of F differ only by a term 
of order 1. Namely, for some C > independent of zq 
and V, 



V 



d^z\{z\zo)f[f{zV-'/')-f{zoV-'/')) 



< C. 



Hence, in particular, p(/i, A,e) = p"(fi,X,e) in the TL. 
Moreover, if /(C) = for IC — CI ^ then the pressure 
is independent of e for ]e| small enough (depending only 
on S). This can be seen as follows. We have 



p"(M,A,e)-p"(/i,A,0) 



1 



■In 



, (16) 



where the last expectation is in the H'^^ ensemble at 
e = 0. The corresponding distribution is exponentially 
localized at z/\/V = C [iTl Il8|. and therefore the right 
side of H16() goes to zero in the TL for small enough 
e. In particular, the e-derivative of the TL pressure at 
e = is zero. By convexity in e, this implies that the 
derivative of p at finite volume, given by (F) f^^x = 
J (Pzf{zV-^^^)Wf,^x{z), goes to zero in the TL. Since / 
was arbitrary, V J^^^Q^g d\ W^,a(Cv^) Q asV ^ co. 



This holds for all 5 > 0, and therefore proves the state- 
ment. 

Our method also applies to the case when the pressure 
is not differentiable in A (which is the case at A = in the 
presence of BEG). In this case, the resulting weights W^,a 
and W^';,^ need not be (5-functions, but Griffiths' method 
\VA Ha] implies that they are, for A ^ 0, supported on 
the real axis between the right and left derivative of p 
and, for A = 0, on a disc (due to the gauge symmetry) 
with radius determined by the right derivative at A = 
0. This, together with convexity, implies that (no) a is 
monotone increasing in A in the TL and, in particular, 

limy^oo V^"^(»^o)m,a=o < limA^o limy^oo ^"^Kao)^',^]^, 
a fact which is intuitively clear but has, to the best of our 
knowledge, not been proved so far in this generality. 
In fact, the only hypothesis entering our analysis, apart 
from the bound on the potential, is the existence of the 
TL of the pressure and the density. 
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